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1 . 


INTRODUCTION 


The progress achieved In the second year of this three-year research pro- 
gram is reported. The ovorall research program is organized into three differ- 
ent subjects: (1) the Jupiter system, (2) the Saturn system, and (3) Comets, 

In the Jupiter system, the untended atmospheres of Xo and their relations and 
interactions with the planetary magnetosphere are the topic of interest. In 
the Saturn system, the hydrogen torus of Titan, the atmospheres associated with 
the planetary riru^s, and the extended atmospheres of the E-ring satellites are 
the topics of concern. For comets, the dust and gas clouds in the outer coma 
are the areas of focus. 

The three primary goals of the overall program are (L) to better under- 
stand the characteristics of the locul atmosphere of lo and the manner in which 
it interacts and delivers a significant source of heavy-ion plasmft to the plan- 
etary magnetosphere, (2) to understand the spatial distribution of hydrogen in 
the Saturn system In terms of thu planet, its rings and its satellite, and to , 
better define the source of protons for the magnetosphere, and (3) to provide 
fresh insight into the chemical and physical nature of comets through study of 
their extended atmospheres. These research objectives are important in achiev- 
ing better understanding of the basic satellite and magnutospheric processes in 
the outer solar system and basic comet and solar radiation processes in the 
inner solar system. To pursue these goals, the research program has employed 
both exploratory calculations to uncover important physical mechanisms and 
more-refined model calculations to evalute the quantitative significance of 
thesu mechanisms. 

For the Jupiter system, the overall three-year research plan is summarized 
in Table 1. There are three major subjects. The first subject, the Io sodium 
cloud, is divided into five subtopics. Efforts during the second year have 
been focused primarily upon studying the 2-D intensity morphology and the 
directional features of the Io sodium cloud using our new model which incorpor- 
ates the oscillation of the Xo plasma torus about the satellite plane. ' 
Progress in these studies is discussed in Section II, Minor efforts have also 
been expended in initiating the development of thu Io potassium cloud model 
using the appropriate information obtained in the first year of this program. 
This development, which is relatively straightforward, is to be completed in 
the third year and will not be discussed further in this report. The third 
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subject of Table 1 has remained inactive during the second year and will be 
activated in the third year in the event that new gas clouds (other than 
sulfur) are discovered. Exploratory modeling for the lo sulfur cloud (recently 
reported as detected by Durrance et al. 1983) is currently being pursued under 
a separate NASA contract. 

For the Saturn system, the overall three year research plan is summarized 
in Table 2 and Is organized into the study of four objects, each known or 
expected to provide a source of hydrogen atoms and ultimately protons to the 
circumplanetary space. During the second year, major emphasis has been placed 
upon better understanding the role of Titan (and ultimately non-Titan sources) 
in the new light that the planet may itself also provide a significant source 
of hydrogen to the Saturn system (Sheraansky and Smith, 1982 ). Progress resul- 
ting from these efforts is discussed in Section III. 

For Comets, research efforts have been mostly inactive in the second year 
and are scheduled to be activated in the third year. Coraetary models developed 
at AER in previous NASA contracts will be refined in the third year of this , 
research program and applied to comet coma observations. 
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II. THE JUl’ITEU SYSTEM: MODELING THE 10 SODIUM CLOU!) 

2.1 Introduction 

The three-year research plan for the To sodium cloud in Table 1 is subdi- 
vided into general model development and four separate subtopics organized 
around the different types of observational data that are available. The 
second year of effort has been built upon the primary model-development activ- 
ities of the first year and has focused upon an initial phase of model applica- 
tion for the subtopics (i) and (iv) of Table 1. The improved model, which now 
includes the oscillating lo plasma torus as a time-space varying sink for 
sodium (see Figure 1), provides a new and powerful modeling resource that has 
already allowed us to draw some significant conclusions from observational data 
for these two subtopics in the second year. These conclusions are discussed 
below. In addition, the computer code for the model has been improved in the 
second year by rearranging Information inside of the CRAY 1 vector loops so 
that the cpu execution speed is now about 30% faster. , 

2.2 The Region B Cloud 

The new sodium cloud model, that explicitly includes the oscillating 
plasma torus, has been used to calculate the two-dimensional intensity morphol- ( 
ogy of the brighter portion of the cloud (i.e., the Region B cloud) on the sky 
plane for several of the sodium cloud images obtained by Murcray (1978). One 
of these observed images at a satellite phase angle of 129° and a system III 
magnetic longitude of lo of 9° (based on the 1965 system in which the north 
pole is inclined toward 200°) is shown in Figure 2. A model calculation for 
these same conditions is shown in Figure 3. 

In the model results in Figure 3, sodium atoms were emitted radially and 

isotropically from Io's exosphere (assumed 2600 km in radius) with an initial 

*•1 ft —2 

speed of 2,6 km sec” 1 and with an assumed constant flux of 4.8 x 10 atoms cm 

sec”* (based upon the surface area of the satellite). The effects of solar 
radiation pressure (although included as an option in the model) were excluded 
in this calculation as a first step in evaluating the impact of the plasma 
torus on the cloud. These effects are expected to be first order corrections 
to the central conclusions to be drawn here and will be included in the explor- 
atory calculations considered in the third year. In Figure 3, the overall 
shape and absolute D 2 Intensity ‘calculated for the cloud beyond the occulting 
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RADIAL DISPLACEMENT FROM JUPITER (Rj) 


Figure li Sodium Electron Impact Ionization Lifetime in Io Plasma 
Torus. The lifetime calculation is based upon the 
measured cross-section by McFarland (1965) , McFarland 
and Kinney (1965) and Zapesochnyi and Aleksakhin (1968) , 
the Voyager 1 data for the plasma density (Bridge, 

Sullivan and Bagenal, 1980; Bagenal and Sullivan, 

1981) and for the electron temperature deduced from 
the Voyager 1 Utfs data (Shemansky, 1980) and jLn situ 
measurements (Scudder, Sittler and Bridge, 1981). 
i The two-dimensional common temperature model of Bagenal 
and Sullivan (1981) , was selected to describe the electron 
density. 
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Figure 2. Io Sodium Cloud Image. The D 2 intensity image was measured in 1976 by Murcray . (1978) 
with a central rectangular mask covering Io and the brightest intensity near the 
satellite. The outer. contour brightness is 0.5 kR (kilo Rayleigh) with contours 
increasing inward by 0.5 kR steps. lo's orbit projected onto the sky plane is 
shown by the dashed line. 
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mask are similar Co Che observed image in Figure 2. Behind Che occulCing mask 
in Figure 3, the cloud intensity increases with decreasing distance from Io by 
almost two orders of magnitude. This sharp gradient in the cloud intensity is 

l 

consistent with the large Integrated D 2 Intensity values measured by Bergstralh 
et al. (1975, 1977) through a 3 x 8 arc sec reactangular aperture centered on 
lo and with unpublished values of Che 1)2 intensity in a 9 arc sec diameter 
circle centered on Io measured more recently by Goldberg (1982). This result 
is particularly gratifying and provides a natural solution to our earlier non- 
oscillating torus modeling results [see: NASA Final Report for Feriod February 

1979 to March 1980, page 46) where the central intensity near Io was much too 
large. The oscillating of the plasma torus, as expected from Figure l, reduces 
the effective strength of the sink more for atom trajectories that populate the 
forward cloud (i.e., atom trajectories inside of Io's orbit and also Inside of 
the ionization sink of the plasma torus) than for atom trajectories that move 
outside and behind Io's orbit. 

Further comparisons of model calculations with other observed Images of , 
Murcray also provide similar results to that of Figures 1 and 2, so that three 
very important conclusions may be drawn: 

(1) The presence (at non-critical Io phase angles) of a predomi- 
nantly forward sodium cloud and the absence of a predomi- 
nantly trailing sodium cloud is consistent with an Isotropic 
emission of sodium from Io and the radial structure of an 
oscillating plasma torus which produces a very asymmetric 
electron-impact ionization sink for sodium (see Figure l), 

(2) The length of the forward sodium cloud is determined natur- 

ally by the plasma torus and occurs when the orbital paths of 
sodium atoms in the forward cloud have their second ioniza- 
tion encounters with the plasma torus somewhat ahead of Io 
(see Figure 4, which also illustrates how the forward cloud 
length is expected to be further altered by the effects of 
splar radiation pressure). ' 

(3) The flux of sodium from Io necessary to produce the measured 
two-dimensional intensity morphology of the B-cloud and also 
the much larger intensity values very near Io is approxira- 

Q rt — 1 

ately 5 x 10° atoms cm 4 sec A , (i.e., a total source of 
about 2 x 10 2 ^ atomd sec -1 , with a more refined value to be 
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^ * Interaction of trie Sodium Cloud and the Plasma Torus. The spatially— projected overlap 
in the satellite plane of the Io plasma torus with the solar radiation perturbed cloud 
shape (solid boundary with heavy shading) and with the unperturbed cloud shape (heavy 
dashed boundary) is indicated for the diametrically opposite satellite phase angles. 

The cloud shape illustrates the envelope of sodium atoms after 20 hours of flight time. 



achieved when a more explicit description of the initial exo- 
spheric velocity dispersion is deduced or otherwise 
measured) . 

i 

Comparisons of additional model calculations with the observational data of 
Murcray (1978) and with more recently acquired but unpublished data of Goldberg 
(1983, with whom a collaborative effort has recently been established) are 
required to address the more refined features of the region B cloud, such as 
the east-west asymmetries (Bergstralh et al. 1975, 1977; Smyth, 1979; Goldberg 
et al. 1980; Goldberg 1981; Smyth, 1983) and the north-south alternating asym- 
metry (Trafton and Macy, 1975; Trafton 1977, 1980), and other small scale 
changes that may appear in the cloud as the result of the oscillation of the lo 
plasma torus about the satellite plane. Such comparisons will provide informa- 
tion about the initial velocity distribution of sodium atoms emitted by lo and 
the nature of the local satellite atmosphere. 

2.3 Interaction with the lo Plasma Torus ( 

Ions and electrons in the lo plasma torus may interact with the sodium 
cloud in two ways. First, the plasma torus determines the lifetime of sodium 
atoms through, the inelastic collisional processes of either electron impact 
ionization (as illustrated in Figure 1) or charge exchange reactions with heavy 
ions. Second, the heavy ions in the plasma torus may provide a non-gravita— 
tional acceleration for some cloud atoms through elastic collisions. Relevant 
inelastic and elastic reactions for sodium in the plasma torus are summarized, 
together with the sodium atom reaction time, in Table 3. The electron-impact- 
ionization reaction time is a minimum value and is taken from Figure 1. Cross 
sections used to compute the charge-exchange reaction times were obtained from 
Johnson and Strobel (1982) and Rutherford et al, (1972), Reaction times for 
the elastic ion-neutral collisioas were taken from Brown, Pilcher and Strobel 
(1983). 

From Table 3, the lifetime of sodium atoms in the lo plasma torus can be 
seen to be dominated by electron impact ionization. Inside of about 5.5 ' 
Jupiter radii, where the electrons become too cool to ionize sodium (see Figure 
1), the dominant plasma torus ions (0 + and S’*") have very small cross sections 
and are incapable of modifying the sodium abundance. Only the effects of elas- 
tic collisions (with reaction time _> 60 hours) and charge exchange with Na^* 
(with reaction time ^ 50 hours) are capable of modifying the sodium abundance 
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Reactions for Sodium Atoms in the Io Plasma Torus 


Sodium Atom 
Reaction Time 

Reaction (hr) 


1. 

e 

+ 

Na 

*► 

Na+ + 2e 


2. 

0+ 

+ 

Na 

► 

0 Na'*" 

W 

3. 

S+ 

+ 

Na 

4 

S + Na + 

00 

4. 

Na + 

+ 

Na 

4 

Na + Na + 

~50 a 

5. 

0++ 

+ 

Na 

4 

0+ + Na + 

~30 b 

6. 

s ++ 


Na 

4 

S+ + Na + 

~20 b 

7, 

hv 

+ 

Na 

4 

Na + + e 

~400 

8. 

X + 

+ 

Na 

4 

X + + Na (elastic collision) 

~20-60 c 


Assumed a nominal ion density of 50 cnf'* (i.e., ~5% abundance) and a 60 km 
sec“* ion-neutral velocity 

Assumed a nominal density of 10 J cm and a 60 km sec 1 ion-neutral ivelo- 

i 

city 

Assumed values between 1000 and 3000 positive charges cm~^ 
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in this region and these reactions are limited vertically to within about 0.1 
to 0.2 Jupiter radii of the centrifugal equator because of the cold ion temper- 
ature. For the region-U cloud, which resides significantly inside of lo’s 

i 

orbit for a time of order 20 hours (see Figure 4), the effects of the plasma 
torus on the sodium cloud abundance inside of 5.5 Jupiter radii should be 
small. The lifetime of sodium determined by electron impact (see Figure 1) is 
therefore the only collislonal process in the plasma torus that is required to 
model the sodium abundance realistically. Ion-neutral collislonal interactions 
of the plasma with the satellite surface and its local atmosphere appear to be 
important, however, as is evident in the discussion of the directional features 
in the next section. 

2.4 The Directional Features 

The new sodium cloud model has also been used to study the directional 
features of the sodium cloud discovered by Pilcher (see Hartline, 1980) in data 
taken at the Mauna Kea Observatory in 1980 and 1981. To understand these f ea- , 
tures, a collaborative effort with Pilcher has been in progress during the last 
two years* During the past year, the basic idea under evaluation has been to 
determine if the directional features, observed outside of lo's orbit on the 
sky plane, could be understood in terms of a high velocity (~20 km sec”*) 
source and the time-varying sodium sink produced (through electron-impact ion- 
ization) by the plasma torus as it oscillates about the satellite plane. 

To investigate this idea a correlation analysis procedure was devised and 
implemented to analyze the relationship between the model-calculated, time- 
varying D-line intensity probability on the sky plane and the vector orienta- 
tion of sodium atoms emitted from Io. This procedure was then applied to study 
several sky plane images of the directional features supplied to us by 
Pilcher. These images showed that the features appearing on the sky plane out- 
side of lo's orbit could change their direction from south to north of the sat- 
ellite plane in a period of several hours. The results of applying our corre- 
lation procedure to these images was to show that to produce a feature thht 
changes its direction from south to north with the observed value of the Io 
magnetic longitude, the atom emission velocities from Io would be required to 
have some limited angular dispersion that is approximately centered at right 
angles to the direction of motion of the satellite. This direction of emission 
velocities is the same as determined by Dr* Robert E. Johnson (University of 
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Virginia) and his graduate student Edward Sieveka (Johnson and Sieveka, 1982) 
who have shown, through elastic collisional calculations for moving ions and 
stationary sodium atoms, that sodium is preferentially scattered at nearly 
right angles to the ion flow direction. 

Model calculations performed near the end of this second program year have 
confirmed the results of our earlier correlational analysis. From this we con- 
clude that the directional features provide substantial support for a satellite 
gas escape mechanism that is driven by the planetary magnetospheric wind as it 
flows around lo. The results of this analysis are far reaching and will be 
reported in a paper to be completed early in the third year of this program. 
These results will also be reported at IAU Colloquium No. 77, "Natural 
Satellites" (Pilcher, Smyth and Combi, 1983). 


/ 


19 



III. 


Tllli SATURN SYSTEM; 


UNDERSTANDING TUB DISTRIBUTION OF 11-ATOMS 


3 . 1 Ini roductlon 

For the Saturn system, a cc^densed summary of our three year research pro- 
gram is given in Table 2, Of the four subtopics listed in Table 2, the Titan 
hydrogen torus has been a long-standing topic at AEK (Smyth, 1981) which was 
supported by NASA several years before this current program began. Efforts at 
AER to understand the extended ring and E-ring satellite atmospheres were ini- 
tiated in 1980-81, one year before the current program began. For this reason, 
efforts in the first year of this research program (see Table 2) were directed 
primarily to the planetary ring utmosphere, where new data had become available 
from the Voyager 1 and Voyager 2 encounters. In the second year of this pro- 
gram, however, the planet was recognized us being a possible source of H for 
the larger Saturn Bystem (Shemansky and Smith, 1982) and the scope of the prob- 
lem has been appropriately broadened und refocused us indicated by the addition 
of the fourth topic in Table 2. An extended hydrogen atmosphere for Saturn 
must now be considered as a plausible source of H atoms previously interpreted 
as a planetary ring atmosphere and might even be an important source for some 
of the hydrogen previously associated with u Titan source. 

At present, therefore, there appear to be four likely sources to supply 
the neutral hydrogen distribution observed in Saturn's magnetosphere. To fur- 
ther explore and understand this subject, our past collaborative effort with D. 
E. Shemansky (regarding the Titan torus data obtained by the UVS instrument of 
the Voyager spacecrafts) has been expanded in scope thiB year. The first step 
in the cxpunded collaboration with Shemansky has been to begin to update the 
Titan model to include the lifetime processes (electron impact ionization, 
charge exchange, and photoionization) for hydrogen operative in Saturn's magne- 
tosphere. Model calculations of the Lyman-a intensity in the magnetosphere 
produced by H atom loss by Titan may then be properly calculated and compared 
to the UVS measurements of the Voyager spacecraft. By studying differences in 
these calculated and measured intensities, the importance of non-Titan gad 
sources may be assessed. Such analysis may allow us to determine, for example, 
if gases from the E-ring satellites Tethys and Dione (for which spatial gas 
envelopes are illustrated in Figure 5) are relevant in understanding the 
Voyager UVS data. Specific efforts in the second year to improve the Titan 
hydrogen torus model are discuss'ed below. 
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Figure 5. Spatial Envelopes for Neutral Gas Clouds of Tethys and Dione. The spatial envelope 
that contains gas atoms emitted from the satellite with velocities of 0.5, 1.0 and . 
1.5 km/sec. is shown in relation to Saturn, the .main planetary rings and other E-ring 
satellites . 


3.2 The Role of Titan 

Efforts in the past year toward understanding the role of Titan in the 
production of Saturn's hydrogen torus have centered around (1) comparing the 
previous modeling work dona at AER (Smyth, 1981) with Voyager 1 and Voyager 2 
data, (2) gathering the information necessary to specify the plasma conditions 
operative at Saturn so the spatial dependence of the hydrogen lifetime can be 
properly modeled, and (3) updating the previous Titan hydrogen model computer 
code to include the spatially varying sink, Progress in this areas has been 
recently reported at the Fifth Conference on the Physics of the Jovian and 
Saturnian Magnetospheres (Combi and Smyth, 1983). 

Figure 6 shows the previously modeled Lyman a distribution on the sky 
plane from a Titan hydrogen source with a flux of 3 x IQ 7 cm sec A from an 
exosphere of 5000 km radius and an ejection velocity of 2 km sec”* (Smyth, 

1981). A constant lifetime of 2 x 10 7 sec had been assumed for hydrogen 
atoms. Based on the results of the Plasma Science experiment on Voyager 2, 
Bridge et al. (1981) calculated the radial dependence in the equatorial plane . 
of the hydrogen lifetime. The calculations Included electron impact ioniza- 
tion, charge exchange, and photoionization processes. Figure 7 compares the 
spatially dependent lifetime with that used in the model. Although scale 
height effects will cause the lifetime to increase with distance from the equa-, 
tor, the major effect of this type of lifetime structure will deplete hydrogen 
originating at Titan inside of 10 Saturn radii (Rg) as compared with the model 
run. 

The principal direct observational data which can be compared with the 
model come from the Ultraviolet Spectrometers (UVS) on board Voyager 1 and Voy- 
ager 2, The Voyager 1 UVS scanned the hydrogen torus along the equatorial 
plane in the light of Lyman a, Broadfoot et al. (1981) concluded that a uni- 
form distribution of hydrogen (10 cm" 3 ) in a cylindrical torus from 8 to 25 Rg 
in radius and at least 6 Rg in height could reproduce the scan data. Figure 8 
shows the UVS scan compared with a similar scan recently produced from the ori- 
ginal model results. Note that in the model 3can the Titan source can be 
directly seen at eastern elongation; also note that there is a definite 
enhancement at Titan's orbital distance (20 Rg) even at western elongation. 

This enhancement is not clearly seen in the UVS data. In the UVS data, Titan's 
location around Saturn is not fixed at eastern elongation, but has an angular 
distribution around the planet that has not yet been determined from the data 
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Radial Displacement from Saturn (Rg) 
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the lifetime calculated By Bridge et al. (1981) from the Voyager 1 PLS data 







DISTANCE FROM SATURN (R s ) 



20 0 20 


DISTANCE FROM SATURN <R S ) 

Figure 8. Model vs. UVS Lyman ot Scan of the Hydrogen Torus. Below 
is the Voyager 1 UVS Lyman a scan, the solid line is the 
simple constant density torus discussed in the text. The 
central peak in the data points is due to the convolution 
of Saturn’s Lyman a disk by the UVS slit. Above is the 
corresponding scan produced from the original model by 
Smyth (1981) for a Titan produced Hydrogen Torus. 
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set. Another difference in Figure 8 between the data and model is the decrease 
in brightness moving in from 20 Rg in the model but the continued increase in 
brightness in the data all the way to 8 Rg. 

A slightly different comparison is shown in Figure 9. This shows the num- 
ber of hydrogen atoms per .1 Rg cylindrical shell as a function of radial dis- 
tance in the equator for both the Titan model and the simple UVS constant den- 
sity model. The underestimation of the hydrogen abundance from 8 to 13 Rg in 
the model can be clearly seen. 

Sandel at al. <1982) reported the Voyager 2 UVS Lyman a scan perpendicular 
to the equatorial plane which showed that the hydrogen torus extends to 6-8 Rg 
above and below the equatorial plane. This would require a Titan hydrogen 
source velocity of 2.6 km s"'* from its exosphere which is also consistent with 
the temperature of 186 K reported by Smith et al. (1982) for Titan's 
exosphere. This higher initial velocity would also populate radial distances 
within 6 Rg from Saturn. 

The question left unanswered at this point is whether the Increased velo-, 
city Implied by the Voyager 2 data will increase the hydrogen density inside 10 
R g enough to offset the effects of a severe sink implied in this region by the 
plasma conditions. If the sink is too strong, it may be necessary to consider 
a source or sources in addition to or instead of Titan for Saturn's hydrogen 
torus. 

In order to answer this question it is necessary to describe quantita- 
tively the spatial dependence of the hydrogen lifetime in two dimensions: the 

radial distance from Saturn in the equator atw the vertical displacement from 
the equator. This effort in collaboration with D. E. Shemansky was discussed 
in the previous section. Some preliminary but as yet incomplete plasma data 
have been obtained from the Voyager PLS team (McNutt, 1983) in this regard. 

Finally, the structure of the previous Titan hydrogen model computer code 
has already been successfully updated to include the two-dimensional lifetime 
structure, once the data are available. At that point only the correct life- 
time data needs to be added in order to completely update the code. 
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